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Nucleation of branes by a four-form field has recently been considered in string 
motivated scenarios for the neutralization of the cosmological constant. An interest- 
ing question in this context is whether the nucleation of stacks of coincident branes 
is possible, and if so, at what rate does it proceed. Feng et al. have suggested 
that, at high ambient de Sitter temperature, the rate may be strongly enhanced, 
due to large degeneracy factors associated with the number of light species living 
on the worldsheet. This might facilitate the quick relaxation from a large effective 
■ cosmological constant down to the observed value. Here, we analyse this possibility 

in some detail. In four dimensions, and after the moduli are stabilized, branes inter- 
act via repulsive long range forces. Because of that, the Coleman-de Luccia (CdL) 
instanton for coincident brane nucleation may not exist, unless there is some short 
range interaction which keeps the branes together. If the CdL instanton exists, we 
find that the degeneracy factor depends only mildly on the ambient de Sitter tem- 
perature, and does not switch off even in the case of tunneling from flat space. This 
would result in catastrophic decay of the present vacuum. If, on the contrary, the 
CdL instanton does not exist, coindident brane nucleation may still proceed through 
a "static" instanton, representing pair creation of critical bubbles - a process some- 
what analogous to thermal activation in flat space. In that case, the branes may 
stick together due to thermal symmetry restoration, and the pair creation rate de- 
pends exponentially on the ambient de Sitter temperature, switching off sharply as 
the temperature approaches zero. Such static instanton may be well suited for the 
$h ' "saltatory" relaxation scenario proposed by Feng et al. 

I. INTRODUCTION 
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It has long been recognized that the effective cosmological constant A e // may have con- 
tributions from a four- form field F, and that in such case 

Kff = (F 2 /2) + A (1) 

may vary in space and time due to brane nucleation events. This has led to various proposals 
for solving the cosmological constant problem, starting with the pioneering work of Brown 
and Teitelboim These authors considered a cosmological scenario where A e jf is initially 
very large and positive, due to a large F 2 term. The additive constant A in (£[]) is assumed 
to be negative, but not fine-tuned in any way, so its absolute value is expected to be of the 



order of some cut-off scale to the fourth power. During the cosmological evolution, A e ff is 
"neutralized" through successive nucleation of closed 2-branes (charged with respect to the 
form field), which decrease the value of F, until eventually A e ff is relaxed down to the small 
observed value, A b s . 

One problem with the original scenario is that neutralization must proceed in very small 
steps, so that any initially large A e ff can be brought to A f, s without overshooting into 
negative values. For that, the charge of the branes should be tiny, ensuring that AA e // < A b s 
at each step. Also, the nucleation rate must be very small, or else the present vacuum 
would quickly decay. These two constraints make the relaxation process extremely slow 
on a cosmological time-scale. Meanwhile, ordinary matter in the universe is exponentially 
diluted by the quasi-de Sitter expansion, resulting in a disappointing empty universe. 

Recently, Feng et al. (FMSW) |2| have suggested that nucleation of coincident branes 
may offer a solution to the "empty universe" problem. Their proposal can be summarized 
as follows. In the context of M-theory, a stack of k coincident D-branes supports a number 
of low energy degrees of freedom, corresponding to a U(k) super Yang-Mills (SYM) theory 
living on the world-sheet. Consequently, the nucleation rate of coincident branes should be 
accompanied by large degeneracy factors, and could in principle be enhanced with respect 
to the nucleation of single branes. The charge of a stack of branes can be very large even 
if the individual charges are small, facilitating quick jumps from A e ff to A t, s . In this way, 
neutralization might proceed very rapidly, perhaps in just a few "multiple" steps of the right 
size. Finally, the stability of the present vacuum could be due to gravitational suppression 
of the nucleation rate [J, |j| . 

FMSW argued, rather heuristically, that the nucleation rate of coincident branes should 
be enhanced by a factor of the form 

D ~ e s , (2) 

where S is the entropy of the worldsheet SYM fields. This entropy was estimated through 
simple thermodynamic arguments, as 

S ~ g*R 2 T 2 , (3) 

where g* is the effective number of worldsheet field degrees of freedom, and R is the size of 
the brane at the time of nucleation. However it remained unclear in [2J which temperature 
T should be used for the worldsheet degrees of freedom. Brane nucleation takes place in an 
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ambient de Sitter (dS) space characterized by a Gibbons-Hawking temperature T Q oc A^^. 
The region inside the closed brane has a smaller value of the effective cosmological constant, 
and is therefore characterized by a smaller temperature Tj. Feng et al. considered two 
alternative possibilities for the temperature of the worldsheet degrees of freedom: T\ ~ T Q 
and T2 ~ (ToTj) 1 / 2 . The proposed enhancement of the nucleation rate and the resulting 
cosmological scenarios are quite different in both cases, and therefore it seems important to 
try and clarify the issue of which temperature is the relevant one. 

The purpose of this paper is to present a more formal derivation of the nucleation rate 
corresponding to multiple brane nucleation. As we shall see, the temperature relevant for the 
worldsheet degrees of freedom is in fact determined by the internal geometry of the world- 
sheet For the Coleman-de Luccia (CdL) instanton, this worldsheet is a 2+1 dimensional 
de Sitter space of radius R, and the corresponding temperature is T ~ R" 1 . When substi- 
tuted into the naive expression (JHJ), this leads to S ~ g*, independent of R (and hence on the 
ambient dS temperatures). As we shall see, the actual result has a certain dependence on 



R due to the anomalous infrared behaviour of light fields in the lower dimensional de Sitter 
space, but this results only in a rather mild dependence on the ambient dS temperatures. 

We shall also see that de Sitter space allows for a "static" instanton which may be 
quite relevant to the nucleation of coincident branes. This is analogous to the instanton for 
thermal activation in flat space. It has a higher Euclidean action than the CdL solution, and 
hence (ignoring the degeneracy factor) it seems to represent a subdominant channel of decay. 
However, we shall argue that, depending on the short distance behaviour of the interactions 
amongst the branes, the CdL instanton for coincident brane nucleation may simply not exist, 
and in this situation the static instanton may be the relevant one. In several respects, the 
static instanton appears to be better suited to the neutralization scenario proposed by Feng 
et al. than the CdL one. 

The paper is organized as follows. In Section II we review different proposals for neutral- 
ization of A e ff via brane nucleation. Section III contains a discussion of coincident branes 
in 4 spacetime dimensions. These are obtained from dimensional reduction of type IIA su- 
pergravity in ten dimensions. In the 4 dimensional picture, the gravitational and four-form 
forces are both repulsive. However, the two are exactly balanced by the attractive force 
mediated by the scalar dilaton. In Section IV and V we discuss the stabilization of the 
dilaton, which is required in a more realistic scenario. After the dilaton aquires a mass, the 
remaining long range forces are repulsive, rendering the stack of coincident branes unstable, 
or metastable at best. This has important implications, since the instanton for nucleation 
of coincident branes will only exist provided that some mechanism causes an attractive 
interbrane force at short distances. 

Section VI contains a description of the CdL instanton for nucleation of coincident branes, 
highlighting a few limiting cases of interest. In Section VII we discuss the corresponding 
degeneracy factor in the nucleation rate, and we show that its dependence on the ambient 
de Sitter temperatures is rather mild. We also include a heuristic interpretation of this 
results based on the observation that the relevant temperature for the worldsheet degrees 
of freedom is determined by the inverse of the radius of the instanton. Implications for the 
scenario of are briefly discussed. 

Section VIII, is devoted to a study of the "static" instanton, where the worldsheet has 
the topology S 2 x S 1 , and where the intrinsic temperature is comparable to T a . In this 
case, the dependence of the nucleation rate on the ambient dS temperature is exponential. 
Coincident brane nucleation can be unsuppressed at large T Q but strongly suppressed at 
present. Our conclusions are summarized in Section IX. Some technical discussions are left 
to the Appendices. 

To conclude this Introduction, a disclaimer may be useful. For most of the paper, we 
shall work directly in four dimensions, and our discussion will be certainly less than rigorous 
from the string theory point of view. In particular, we shall model the degrees of freedom of 
a stack of k coincident 2-branes by a weakly coupled U(k) gauge theory on the world-sheet. 
This may or may not correspond to a true dimensional reduction from M-theory, but it 
should at least represent some of the broad features of the degeneracy factors. 

II. NEUTRALIZATION v.s. RANDOMIZATION 

In four dimensions, a four form can be written as T^pa = F^—gep Vpu where g is the 
determinant of the metric and e is the Levi-Civita symbol. This field has no propagating 
degrees of freedom, since in the absence of sources the equation of motion d*T = implies 



that F is a constant. This simply gives a contribution to the effective cosmological constant 
which gets added to the true cosmological constant, or vacuum energy density A, 
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Ae// = -y + A- (4) 

Four-forms may couple to "charged" 2-dimensional extended sources, or 2-branes, through 
a term of the form 

q J A, (5) 

where the integral is over the worldsheet of the extended object and A is the 3-form potential 
{T = dA). In this case, F changes by 

AF = q (6) 

accross a brane of charge q. Consequently, F can decay through nucleation of closed spherical 
branes. The process is analogous to pair creation in the presence of an electric field, and 
very similar to false vacuum decay in field theory p] . The closed brane is the boundary of a 
newly formed "true" vacuum bubble, where the field strength differs from the original value 
by the amount (JHl). After nucleation, the radius of the bubble grows with constant proper 
acceleration, and the volume occupied by the new phase keeps increasing. Further nucleation 
events take place in the region with a low F, lowering A e ff even further. In the absence of 
gravity, this "neutralization" would proceed as long as F > q/2, wiping out any large initial 
value of F 2 and leaving us with a large negative cosmological constant bJ^ 1 ~ A + 0(q 2 ). 
Of course, this is not what we want. 

Gravitational effects can improve the picture dramatically 0]. In particular, tunneling 
from flat or Anti-de Sitter space (with a negative cosmological constant) is forbidden pro- 
vided that the squared tension of the branes is sufficiently large compared with the jump 
AA e jj ~ qF in energy density across the brane, 

a 2 > (4/3)qFM 2 . (7) 

Here a is the brane tension and M 2 = 1/ (87rG) is the square of the reduced Planck mass. 
This gravitational "shutdown" of brane nucleation could be useful, since an initially large 
K e ff may eventually get "neutralized" to a value which is much smaller than A in absolute 
value Suppose that the true cosmological constant is negative A < 0. As long as 
A e ff > 0, branes keep nucleating. But once a vacuum with A e // < is reached, the process 
stops provided that Eq. (|7J) is satisfied. After that, the vacuum becomes absolutely stable. 
Brown and Teitleboim conjectured that we may live in one such vacuum, where the effective 
cosmological constant is expected to be of the order of the energy density gap between 
neighboring vacua 

\A f jf\ ~FAF~g|A| 1/2 . (8) 

In this vacuum, \F\ ~ (A) 1 / 2 , and the huge negative bare cosmological constant is almost 
completely cancelled by the F 2 contribution in the final state. 

For sufficiently small charge q, Eq. (jHJ) leads to a suppression of |A{^ a/ | relative to the 
true vacuum energy |A|, which might be helpful in solving the "old" fine tuning problem of 
the cosmological constant. Particle physics models suggest |A| > (TeV) 4 . Hence, we need 



q < 10- 35 (eV) 2 , 



(9) 



so that the final value |A^ aZ | is consistent with the observed value |A°jj| ~ 10~ u (eV) A . The 
constraint (jHJ) seems rather demanding, since in the context of supergravity we would expect 
q to be closer to the Planck scale. This is the so-called "gap problem". FMSW argued that 
the smallness of the charge could be due to the wrapping of branes on degenerating cycles in 
the extra dimensions. A successful implementation of this idea has not yet been presented, 
but some plausibility arguments have been given in Alternatively, in a different context, 
it has been suggested that branes with a very tiny charge q may arise due to symmetries 
of the theory. An explicit example was given in [6j], where the branes are not fundamental 
objects but domain walls of a broken discrete symmetry. This same symmetry suppresses 
the coupling of the domain walls to the four-form F without any fine-tuning of parameters 
(see also for a fuller discussion). 

A more severe problem of the Brown and Teitelboim neutralization scenario, is the "empty 
universe" problem which we discussed in the Introduction. By combining the condition 
(JHJ) with the stability condition (JZJ), it can be shown that the time required to reach the 
value (jHJ is huge compared with the age of the universe 0. By the time the effective 
cosmo logical constant would be wiped out, all other forms of matter would have also been 
diluted exponentially, in clear contradiction with observations. Furthermore, the endpoint 
of neutralization would be a state with vanishing or negative effective cosmo logical constant, 
whereas the observed value A°^ is positive. 

One way around the empty universe problem is to consider a slightly different scenario, 
where the effective cosmological constant is "randomized" (rather than neutralized) during 
inflation. Assume, for simplicity, that the energy scale of inflation is much larger than 
AA e ff. In an inflationary phase, brane nucleation processes may increase as well as decrease 
the value of F 2 jj] . Thus, A e jf will randomly fluctuate up and down the ladder as a result 
of brane nucleation. Inflation is generically eternal to the future and there is an unlimited 
amount of time available for the randomization process to take place before thermalization 
jij. Assume that the tunneling barriers are sufficiently high, so that no nucleation events 
happen in the last 60 e-foldings of inflation, or during the hot phase after thermalization 
up to the present time. In this scenario there is no empty universe problem: the local value 
of A e ff is decided many e-foldings before the end of inflation, and a wide range of values of 
A e ff will be found in distant regions of the universe, separated from each other by distances 
much larger than the present Hubble radius. Some of these regions will just happen to have 
a very tiny A e ff. In combination with anthropic selection effects, this approach may be used 
to explain the smallness of the observed effective cosmological constant 0, El, H, E| ■ This 
may also explain the so-called cosmic time coincidence, or why do we happen to live at the 
time when an effective cosmological constant starts dominating 0, 0, E| ■ 

Bousso and Polchinski 0] have proposed a somewhat related "randomization" scenario 
which, moreover, does not rely on branes with the exceedingly small charge q satisfying 
In the context of string theory one may expect not just one but many different four- form 
fluxes Fi coupled to branes with different charges Each one of these fluxes is quantized 
in units of the charge, so that F^ = niqi. In this case, the condition for a generic negative 
cosmological constant to be compensated for by the fluxes is \A e ff \ = \ J2i=i( ( li n 'i)/2 + M ~ 
A ofes . The larger the number of different fluxes, the denser is the discretuum of possible 
values of A e ff, and the easier it is to find a set of values of rii such that the above inequality 
is satisfied. A sufficiently dense discretuum is typically obtained provided that the number 
of fluxes is sufficiently large J ~ 100, even if the individual charges qi are Planckian (a 
smaller number of fluxes J > 6 may be enough in a scenario with large extra dimensions, 



where the charges are suppressed with respect to the Planck scale by a large internal 
volume effect). In the Bousso-Polchinski model, A e ff is typically very large, and drives 
an exponential expansion at a very high energy scale. Suppose we start with a single 
exponentially expanding domain characterized by a set of integers {rii}. Whenever a brane 
of type j nucleates in this region, the integer rij will change by one unit inside of the brane. 
The newly formed region will itself expand exponentially, creating a huge new domain. The 
nucleation of further branes within this region will cause an endless random walk of the 
values of {n^}, which will sample the whole discretuum of values of A e jj. Eventually, a 
bubble will nucleate where A e ff is comparable to the observed value. The nucleation of this 
last bubble is still a high energy process, which kicks some inflaton field off its minimum, and 
starts a short period of inflation within this last bubble. This period of "ordinary" inflation is 
necessary in order to produce the entropy we observe, thereby avoiding the empty universe 
problem. Of course, some anthropic input is still necessary in this approach in order to 
explain why, out of the discretuum of possibilities, we live in a vacuum with small A e jf. 

Nucleation of coincident branes would drastically modify the neutralization scenario of 
Brown and Teitelboim, as well as the randomization scenarios sketched above. As proposed 
in 0, in the case of neutralization, this modification may lead to a solution of the empty 
universe problem. In the randomization scenarios there is no such problem, and it is unclear 
whether an enhancement of the multiple brane nucleation rate is desirable at all. This 
enhancement would trigger large jumps in the effective cosmological constant, making the 
calculation of its spatial distribution more complicated than for single brane nucleation. 
Thus, it is of interest to understand the conditions under which multiple brane nucleation 
is allowed, and what are the degeneracy factors which might enhance their nucleation rate 
relative to the nucleation of single branes. Before addressing this issue, it will be convenient 
to present a short discussion of coincident branes. 



III. COINCIDENT BRANES IN 4D 

In the context of string theory, one may consider stacks of k coincident D-p-branes. Each 
brane in the stack has charge q with respect to the form field A. Thus, in four noncompact 
dimensions, a pair of parallel 2-branes repel each other with a constant force per unit area 
given by 

f q = q 2 /2, (10) 

due to the four-form field interaction. In the ten dimensional theory, the repulsive force due 
to A is balanced with other contributions from the closed string sector, such as the graviton 
and dilaton. As a result, there are no net forces amongst the different branes on the stack. 

It is in principle possible to maintain this delicate balance by suitable compactification 
from 10 to 4 dimensions. In 4D, the branes look like domain walls, and their interaction 
through the ordinary graviton leads to a mutual repulsive force given by 

U = 3a 2 /4M 2 . (11) 

Hence, both forces given by (jl()j) and (fTTJ) tend to push the branes apart from each other. 
On the other hand, some of the higher dimensional gravitational degrees of freedom are 
represented by scalars in 4D, and these, together with the dilaton, lead to attractive forces. 



A. Dimensional reduction 



Let us first consider the case of D-2-branes in 10 dimensional Type IIA supergravity. The 
relevant part of the action is given by (see e.g. 
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e~ 2 * (i? + 4(V0) 2 ) 
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T 2 



T 2 / d 3 C V^s e-* + q 2 / d 3 £ A. 
Jt, Jt, 
_ (12) 

Here, \[G and \/Gy, are the determinants of the 10 dimensional metric Gab arid of the metric 

induced on the worldsheet £, respectively, whereas R is the Ricci scalar corresponding to 

Gab- The hats on the four-form, the gauge potential, and the corresponding charge, are 

introduced in order to distinguish them from the four dimensional ones which will be used 

below, and which differ from those by constant normalization factors. Compactifying on a 

Calabi-Yau manifold K through the ansatz 



ds 2 1Q 



e^-^g^dx^dx" + e 2 ^dK 2 , 



where the Greek indices run from to 3, we readily obtain the following four- dimensional 
action: 

Ml , , , _ „ „ , ., .., 
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(13) 



Here, M 2 = Mf Ve, where V§ is the coordinate volume of the manifold K, TZ is the Ricci 
scalar for the metric g, and we have introduced two linear combinations of the internal 
volume modulus ip and the dilaton <f> 



= 2(p- 9i/j, 



a 



0- ip. 



The field a decouples from the branes, and shall be ignored in what follows. In Ref. [2], 
a different expression was given for the dimensionally reduced action, because no modulus 
was introduced for the size of the internal space. However, as we shall see, such modulus is 
necessary for the cancellation of forces among the branes. 

A more direct route to the 4-dimensional theory (|T3|) . which will be slightly more con- 
venient for our discussion, is to start directly from 11-dimensional supergravity and com- 
pactifying on a 7 dimensional internal space. The action in eleven dimensions is given by 
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T 2 



-T 2 / d^VGx + fa / A. 



where T 2 = q 2 = 2irMf 1 . Introducing the ansatz 

ds 2 u = e-^g^dx^dx" + e 2 ^dtt 2 , 
with a Ricci flat internal manifold, we find 



Sa 



M 2 
p 



(14) 



(15) 



(16) 



where M 2 = M^Vr and 



2 v 



Not surprisingly, this has the same form as JEJ), since after all the 10-dimensional Lagrangian 
can be obtained from the 11-dimensional one by compactifying on a circle. 

Eq. (|TH|) is a particular case of the slightly more general action in four dimensions: 

S * = \f rf V£ [^1 [K - (d V ) 2 ] - i e" 2 ^ 2 } -a J d 3 eVT e av + g ^ A. (17) 

The parameter a characterizes the scalar charge of the brane. As we shall see below, lin- 
earizing in ip around tp = 0, it can be easily shown that the scalar force is given by 



f e = -e 2 /2, 
where we have introduced the scalar charge 

e = aa/M p . 



(18) 



(19) 



Thus, from (jllj) and pSjl . the branes will be in indifferent equilibrium provided that 
the following relation holds: 



e 2 — q 2 



3a 2 
4M 2 



0. 



(20) 



From (|19|h this condition can be rewritten as 
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(21) 



Note that the case of Type IIA supergravity discussed above, corresponds to a = a/7/2, 
q = (>/2/Mp) and cr = T 2 [note that ^4 = (x/2/M p )^4]. Therefore (f2"Tj) is satisfied provided 
that T 2 = g 2 , the usual BPS condition. 



B. Multiple brane solutions 

Multiple brane solutions to (fTTj) can easily be constructed provided that (f2"Tj) is satisfied. 
In the bulk, the equation of motion for A leads to 

J 7 , otf3y8 = F e 2av ^fg e Q/ 3 7 5, 

where F is a constant and e is the Levi-Civita symbol. At the branes, this constant jumps 
by the amount 

AF = q 

It can be checked that the remaining equations of motion for the scalar and the gravitational 
field follow from the action 

S * = E Jf x ^ (<^) 2 ] - VM } - ff jf d3 ZVl 6°*, (22) 



where 

Viip) = S e 2 ^, (23) 
and the sum is over the regions with different values of F. With the metric ansatz 

ds 2 = w 2 (z)r] ab dx a dx b + dz 2 , 

where Latin indices run from to 2, the solution is given by 

and 

w(z) = e~2«. (25) 

Here, q are integration constants. These, and the sign option in (jUJ), must be chosen so 
that the junction conditions for the gravitational field and for the scalar field are satisfied 
at the branes. For the gravity part, the condition is (l4| . 

[K ah ] = -AirGa e a ^ ab , (26) 

where [K ab ] is the difference of extrinsic curvature on the two sides of the brane and is 
the world-sheet metric. In the present case, this reduces to 
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w 



2M p 2 



e av . (27) 



For the scalar field, the junction condition which follows from the equation of motion at the 
brane reduces to 

¥\ = ^ e " (28) 
p 

which is consistent with (12 7|) and (|25p. For instance, a solution with two branes separated 
by a distance d is given by flat space (w = l,<p = 0) in the region between the branes, and 
by 

e~^ = l-a 2 ^- (\z\-d/2) (\z\>d/2) (29) 
p 

in the exterior region (see Fig QJ. The solution corresponds to two branes of charge q 
interpolating between regions with F = q and F = —q, separated by a region with F = 0. 

Note that the solution ()29|) contains a singularity at a distance z ~ M 2 / o?a (where the 
warp factor w and the volume of the internal space vanish). In what follows, however, we 
shall not be interested in flat, infinite branes, such as the ones discussed above. Rather, we 
shall be interested in compact instanton solutions of finite size, in a theory with a stabilized 
dilaton. In this case, the singularities due to linear potentials of the form (|29|) should not 
arise, but the solutions discussed above should remain a good description in the vicinity of 
the branes. Another consequence of stabilizing the dilaton is that the perfect balance of 
forces amongst the branes will be spoiled at distances larger than the inverse mass of the 
moduli, as we now discuss. 



FIG. 1: Configuration of two parallel branes 



IV. STABILIZING THE MODULI 

As it stands, the dimensionally reduced supergravity lagrangian (fP3*|) [or its generaliza- 
tion (jI7j) ] is not useful for discussing the neutralization of the cosmological constant. The 
lagrangian does not include the bare cosmological term A, which is precisely the subject of 
our interest, and the term proportional to F 2 does not behave as an effective cosmological 
constant, but rather as an exponential potential for a modulus ip (which is not flat enough to 
mimick the vacuum energy 1 ). A more realistic model is obtained by introducing a stabiliza- 
tion mechanism which fixes the expectation value of (p, and gives it some mass m. Once <p is 
stabilized the F 2 term does behave as a contribution to A e //. Stabilization is also desirable 
because the dilaton and the radion moduli (corresponding to the size of extra dimensions) 
mediate scalar interactions of gravitational strength, which are severely constrained by ob- 
servations. The study of mechanisms for stabilization is currently an active topic of research 
(see e.g. [l6L Il7| and references therein). Although the details of stabilization will not be 
too important in our subsequent discussion, it may be nevertheless illustrative to have in 
mind a specific toy mechanism (for which we do not claim any rigorous justification in the 
context of string theory). 

The general problem is that the potential ()23)1 has no minimum and leads to a runaway 
dilaton. In order to create a minimum, let us consider two contributions which may be 
added to IJ2Hj) . First of all, instead of using a Ricci flat internal manifold in (|15j) we may 
compactify on an Einstein manifold, with 

Rab = (30) 



Exponential potentials such as the one appearing in (|22fl have been thoroughly studied in the literature, 
and it is known that they can drive cosmological solutions with a power-law scale factor 0|. Such 
attractor solutions are approached for a wide range of initial conditions, and the resulting expansion 
can be accelerating or decelerating, depending on whether a 2 < 1/2 or a 2 > 1/2. For a 2 < 3/2 the 
cosmological scale factor approaches a(t) ~ i 1 / 2 " 2 , where t is cosmological time. This solution corresponds 
to an effective equation of state p = [(4a 2 /3) — l]p, where the ratio of kinetic and potential energies of 
the scalar field ip remains constant. In our case, from 1)21(1. we need a 2 > 3/2 and therefore the kinetic 
term becomes completely dominant in the long run, which leads to p = +p. Hence, by itself, the F 2 term 
does not behave like an effective cosmological constant. 



Upon dimensional reduction, the curvature K contributes an exponential term to the ef- 
fective potential for ip. However, this will still not be sufficient for stabilizing the internal 
volume in an interesting way. In fact, as shown by Nunez and Maldacena [18j . there are 
no static compactifications of the classical supergravity Lagrangian with a positive effec- 
tive four-dimensional cosmological constant A e ff > 0, and so far we have added nothing 
to the classical Lagrangian. Thus, in order to implement the four dimensional situation of 
our interest, a third term related to quantum corrections has to be considered. Following 
Candelas and Weinberg we may consider the Casimir energy of bulk fields (we are of 
course assuming that supersymmetry is broken, so that the Casimir contributions of bosons 
and fermions do not cancel each other exactly). 

In the example considered by Candelas and Weinberg besides the Casimir energy 
term, a higher dimensional cosmological term Ai +n was used, and the internal manifold 
was taken to be a space of constant positive curvature (K = 1). In 11-D supergravity a 
cosmological constant is not allowed, so here we use the F 2 flux instead. Also, we will need 
to compactify on a negatively curved internal manifold (K = —1), or else the dilaton would 
be stabilized at negative A e //. 

The versatility of negatively curved compactifications has been stressed in (20l |. In par- 
ticular, they have the interesting property of rigidity, which means they do not lead to other 
moduli besides the size of the internal space. Compact hyperbolic manifolds (CHM) can be 
obtained from the maximally symmetric negatively curved space H 7 through identifications 
by a discrete isometry group T. The volume of H 7 /Y is given by 

V 7 = rje 7 , (31) 

where r c is the curvature radius of the manifold, related to the curvature parameter in 
(J30J) by K = —1/r 2 . The factor e 7 depends on the topology, and it is bounded below but 
not above. If L is the largest distance around the manifold, then for L ^> r c /2 we have 
e 7 ~ e 6L//rc . The Kaluza Klein (KK) spectrum in this manifold is believed to have a mass 
gap, bounded below by ttlkk ~ e~^r~ 1 . From the 11 dimensional point of view the Casimir 
energy density scales like 

pg 1} = Cm% K = Ce-^\ (32) 

The factor C can be estimated by naive dimensional analysis as C ~ f3v, where v ~ 10 3 is 
the number of physical polarizations of bulk fields, and (3 is some small one-loop factor. This 
factor will depend on the precise topology of the compact hyperbolic manifold, but it could 
plausibly be in the range (3 ~ 10~ 2 — 10~ 4 . Hence, the parameter C could be of order one. 
Explicit calculations for different choices of the manifold H 7 /Y have not been performed, 
and are well beyond the scope of the present paper (see e.g. |2j| and references therein). 
In what follows we shall leave C unspecified, assuming that a compactification exists where 
\C\ > 1. Multiplying the higher dimensional energy density p@ given in ()32j) by the size of 
the internal manifold, V 7 e 7 ^, and by a factor e -14 ^ which arises from the four- dimensional 
volume element in going to the Einstein frame, the effective potential which appears in Eq. 
(J22J gets replaced by 

Vi(<p) = -21 Ml x V 7 Ke~^ + C V 7 r; u e - 18 ^ + ^e" 21 ^, (33) 



where we have also added the curvature contribution. Here ip = — (2/2l)cap, where a = 



We can always adopt the convention that r c = 1/Mn, since a change of r c — > e x r c can 
be reabsorbed by a shift in ip — > ip + A and a constant re-scaling of the four- dimensional 
metric g^ v — > e 7X g fJ-l/ (in this frame, the curvature of the manifold is of the order of the 
higher dimensional Planck scale for \ip\ < 1). Since Fi = n^g, where is an integer and 
q = (\/2/M p )q 2 = 2y/2irM% 1 /M p , we have 



Vi(<p) = M 2 p M 2 u 



21 e-^ + C e 



- 18 ^ + 4ttV 2 



M 



11 



-21-0 



(34) 



As illustrated by Eq. (pTTJl a large value of Vj in units of r c = 1/Mn can be obtained by using 
a manifold with sufficiently complicated topology. Since M 2 = M^V? = M^e 7 , the factor 
(M n /Mp) 4 = e" 27 in the last term of flHU) can be rather small. The scale M n could be as 
low as the TeV, in which case that factor can be as low as 10 . Moreover, as emphasized 
in , this can be achieved for CHM even if the linear size L of the internal manifold is not 
very much larger than r c . In Fig. 2 we plot the effective potential for C = —20, e -7 = 10~ 3 
and values of rii = 485, 487, 489 and 491. 

When a single brane nucleates, it changes rii by one unit, and hence changes the value 
of the effective potential at the minimum (changing therefore the effective cosmological 
constant). If the discretuum of values of Fi were sufficiently dense, then there would always 
be one of the minima of the effective potential where the vacuum energy is sufficiently 
small to match observations. In the case we have considered here, the discretuum is not 
dense enough. The cancelation between the last term in (}3*3*j) and the other two requires 
rii ~ (M p /Mu) 2 , and so the gap between levels near Vi = can be estimated as AVi ~ Mf x , 
which is far too large. This situation can be remedied by considering 5M branes wrapped 
around 3-cycles in the internal space. As emphasized by Bousso and Polchinski [ll[, these 
are coupled (with different charges) to additional fluxes, and a large number of fluxes will 
result in a much denser spectrum. Alternatively, FMSW have suggested that the branes 
may wrap a degenerating cycle Q, in which case the individual charges might themselves 
be exponentially smaller than the fundamental scale, resulting also in a sufficiently dense 
discretuum. 
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FIG. 2: The effective potential (]34[). for different values of the integer m which characterizes the 
quantized flux of the four-form F. 



Before closing this Section, we should note that the location where the modulus sits is 
basically determined by the competition between the two first terms in ()33|) . corresponding to 



curvature and Casimir energy. The physical curvature radius of the internal space is therefore 
stabilized at r p h ys = r c e^ ~ M xx C l ^ D ~ 2 \ where D = 11 is the spacetime dimension. 
Thus, unless the constant C in Eq. (J3~2"j) is exceedingly large, the compactification scale 
is comparable to the inverse of the higher dimensional cut-off scale M n . In such case the 
semiclassical analysis which we have entertained above is not justified, since higher order 
corrections will be just as important as the one loop effect which we have included. This 
appears to be a generic problem when we try to stabilize by making the curvature and 
the Casimir terms comparable, as in the Candelas and Weinberg example 0]. There, the 
problem was not quite as poignant, since the constant C could be made very large by adding 
a sufficient number of fields (also, as it is clear from the above estimate of r p h ys , the problem 
is somewhat milder if the number of extra dimensions is smaller). Here we shall not dwell 
on this problem, since the main purpose of the above discussion is just to illustrate the 
role of the four-form in obtaining a discretuum of states. For more fundamental approaches 
to moduli stabilization in the present context, the reader is referred to [16| and references 
therein. Our ensuing discussion will be largely independent of the details of the stabilization 
mechanism. 



V. INTER-BRANE FORCES 

Inevitably, the stabilization of (p spoils the perfect balance of forces. At distances larger 
than the inverse of the mass of the modulus m -1 , the remaining gravitational and four- form 
interactions are both repulsive, and lead to a linear potential per unit surface of the form 



V(d) 



I I I I ! » ill ' 

2 4M 2 _ 



q l 3a 



d (d > m- 1 ) (35) 



where d is the distance between the branes. 

To investigate the behaviour of the interaction potential at distances shorter than m _1 , 
let us first consider the situation where gravity and the 3-form gauge potential A are ignored, 
and the branes interact only through a scalar field <p of mass m. The action is given by: 

M 2 r r 

S = — ^ J d 4 x [{dtp) 2 + m 2 <p 2 } - a J d 3 ^ e^. (36) 

The solution with a single brane on the plane z = has the cusp profile 

¥> = ¥>„e- mW , (37) 

where (p is a solution of 

e -a V0 = e_ 

^° 2mM 2 2mM p 1 ; 

The energy per unit area of this configuration is given by 

ill 2 



a x = ae aip0 



M z r 

+ -± dz( V ' 2 + m V) = ae aLp0 + mM 2 (p 2 Q . (39) 
For small charge and tension, we have 



e 2 

CTi « a — - — , (e 2 <C 2ma). 
4m 



Due to the scalar field dressing, the effective tension of the brane, denoted by <7i, is smaller 
than the parameter o which appears in the action. This effect becomes more dramatic if we 
place a large number k of branes on top of each other. Since both a and e scale like k, the 
effective tension of the stack is given by 



a k ^k 



a — k-^— 

Am 



(k e 2 < 2ma), (40) 



which grows with k but less than linearly. In the limit of very large k, using (jHEJ) with a 
replaced by ka we have 

^ In 2 (£*L) , (k e 2 » 2ma) (41) 

and so the tension almost saturates, growing only logarithmically with k. This last expression 
should not be taken too literally, however, since in ()36|) we are neglecting gravitational 
effects. As shown above, nonlinear effects of gravity become important at a distance given 
by M 2 /a 2 a k , which in the limit given by (|41|) is smaller than m~ l . Hence the cusp solution 
(jnZj), which has typical width ~ m~ x , will receive sizable gravitational corrections in the 
limit of large k. Nevertheless, it still seems likely that the effective tension of the stack of 
branes will grow with k much slower than linearly. Physically, the reason is that the scalar 
charge of the stack increases with k. According to (|2*Hj) . this means that the cusp in the 
field if on the branes grows stronger, which in turn means that the value of the field </9 on 
the branes has to be further displaced into negative values. Hence, the brane contribution 
to the effective tension kae ar ^° shows only a very modest growth with k, and the tension 
for large k is in fact dominated by the potential and gradient energy of the scalar near the 
brane. 

Let us now look at the interaction potential between two branes separated from each 
other. For simplicity, we shall restrict attention to the case of small scalar charge, e 2 -C 2ma. 
Then, the last term in ()36|) is well approximated by its linearized expression: 



m 2 c r 

S = --JL J d*x [{dip) 2 + mV] - J d 3 ^ {a + eM p <p). (42) 

Placing the two branes at z — ±d/2, the solution for the scalar field has the "Golden Gate" 
profile shown in Fig. 3 

if = %- e- md/2 cosh mz, \z\ < d/2 (43) 

Y mM p 'ii/ \ J 



up = ^- cosh(mci/2) e" m|z| . 1*1 > d/2 

Y mM p K 1 1 ii/ 



The energy per unit area of this configuration, as a function of the interbrane distance, is 
given by 



a 2 (d) = 2a, - ^—e- md . (44) 
2m 

From this expression, and adding the long range contributions from gravity and the four- 
form, the interaction potential per unit surface is given by 



V(d) 



<r 3<r 2 1 , e 2 _,„, 



2 AM?, 



2m 



At short distances, this takes the form 



V{d) = -— + Q 2 d-—-d 2 + ... (45) 
w 2m 4 v ' 



where the parameter 



4M p 2 



(46) 



was introduced in ([20)1 . As shown in Section III, dimensional reduction from the supergravity 
lagrangian with BPS charges gives Q 2 = 0. In this case, the linear term in (|43j) disappears. 
The quadratic term is negative, which means that the stack of coincident branes is un- 
stable and tends to dissolve. Thus the stabilization of the modulus (p seems to make the 
superposition of branes an unstable configuration. 

As we shall see in more detail in Section VII, stacks of branes in marginal or unstable 
equilibrium will not be appropriate for constructing instantons, since in particular these 
would have too many negative modes. The instanton is only meaningful if the branes 
attract each other. The above analysis shows that at the classical level, the branes with 
BPS values of the charges would not attract each other, and consequently the nucleation of 
coincident branes is not allowed at least in the semiclassical description. 

There may be several escape routes to this conclusion. For instance, after supersymmetry 
breaking, the charges of the branes get renormalized, and it is possible that the corrected 
charges satisfy Q 2 > 0. In this case the branes in the stack would attract each other 
with a linear potential. Other mechanisms by which nearby branes attract each other are 
conceivable, but here we shall not try to pursue their study. It should be emphasized, 
however, that this remains an important open question which needs to be addressed in order 
to justify the semiclassical description of coincident brane nucleation. In section VII, where 
we discuss the degeneracy factor, we shall simply postulate that an attractive interaction 
exists at short distances. 




FIG. 3: Profile of a massive dilaton in the presence of two branes. 



VI. COLEMAN-DE LUCCIA (CDL) INSTANTONS 

The brane nucleation rate per unit volume is given by an expression of the form |j| 

T = De~ B , (47) 



where B = Se(I)—Se(B). Here Se{I) is the Euclidean action of the instanton corresponding 
to the decay of the four-form field, and Se(B) is the action for the background solution before 
nucleation. The prefactor D will be discussed in the next section. The formal expression 
(JUj) can be used both in flat space and in curved space. Also, it can be used at zero or at 
finite temperature. The difference is in the type of instanton and background solutions to 
be used in each case. 

In this section, we shall concentrate in the maximally symmetric instanton. In flat space, 
this represents the decay of a metastable vacuum at zero temperature. The Euclidean 
solution can be described as follows At infinity, the field strength takes the value F Q , 
which plays the role of a false vacuum. Near the origin, we have F = Fi = F — q. This plays 
the role of a true vacuum phase. Both phases are separated by the Euclidean worldsheet of 
the brane, which is a three- sphere of radius 

R = ^. (48) 

Here, a is the tension of the brane and 

z = \{F 2 -F?) = q [F -{q/2)} (49) 

is the jump in the energy density accross the brane. The difference in Euclidean actions 
between the instanton and the background solution is given in this case by |H 

B^ » ^fl- ( 5 °) 

If instead of considering a single brane, we are looking at the nucleation of k coincident 
branes, an analogous solution of the Euclidean equations of motion should be considered. 
The only difference is that the effective charge is a factor of k higher, and the effective 
tension is also higher by approximately the same factor. Using (|4()jl . we should replace 

e ^ e h = kq[F - {kq/2)], 



a -> cr fc ps k[cr - (ke 2 /4m)}, (51) 

in Eq. (j5Up for the "bounce" action. As shown in Section V, the approximate form of 
is valid when the second term in the r.h.s. of (|51|) is small compared with the first. In the 
model considered in Section IV, the modulus is stabilized with a mass of order m ~ M n . 
This may be regarded as a conservative upper bound for the general case. Since e = aa/M p , 
the condition (ke 2 /Am) <C o requires 



a 2 \ a J \M 1 

For larger k, gravitational corrections to the brane profile become important on a lengthscale 
comparable to the inverse mass of the dilaton, as we discussed in the previous section. In 
this case the instanton solution will depend on the detailed dynamics of the dilaton near 
the brane. The investigation of these dilatonic instantons is per se an interesting problem, 
which we leave for further research. Here, for simplicity, we shall assume a scenario where 



either a <C Mf 1} due perhaps to the wrapping of branes on a degenerating cycle, or where 
the extra dimensions are relatively large, so that M p 3> M n , and we shall restrict attention 
to instantons where the number of branes is bounded by (|52|). This gives 

B {fiat) _ 27tt 2 k a 4 
fc ~ 2 g»[F -(*g/2)]»- ^ 

Nucleation in flat space is impossible when the number of branes is too large, since otherwise 
we would be jumping to vacua with a higher energy density. Therefore, we must restrict 
to k < 2F a /q. In fact, the minimum value of Fi is achieved for the largest k satisfying 
k < (F /q) + (1/2). Note that the action increases faster than linearly as we increase the 
number of branes k, but we still have 



B Ulat) ^ k B (flat) 

throughout this range. 

When gravity is taken into account, the maximally symmetric instanton was given by 
Coleman and De Luccia 0. It is constructed by gluing two different de Sitter solutions 
(that is, two four-spheres) at the worldsheet of the brane, which is still a three-sphere. The 
instanton is sketched in Fig. 0] The four-spheres have the radii if" 1 and H~ x , where 

rr2 _ A eff ^2 _ A eff ~ 6 ( , 

H ° ~ 3M|' H <~ ^Mp (54) 

are the Hubble rates of the de Sitter phases before and after the nucleation event, respec- 
tively. Here, and in what follows, we are assuming that the effective cosmological constant 
is still positive after nucleation, since these are the final states we are interested in. The 
bubble radius at the time of nucleation (which coincides with the radius of the three-sphere) 
is bounded by < R < H~ x . Analytic expressions for R are given in Refs.0,0|. The general 
expression is cumbersome and not particularly illuminating, so we shall concentrate on a 
few limiting cases of interest. 




FIG. 4: Coleman-de Luccia instanton, which is obtained by gluing together two different four- 
spheres of radii H^ 1 and H~ l at the worldsheet of the brane, itself a three-sphere. 

As discussed in Section III, the gravitational field of a brane is repulsive, and is character- 
ized by an "acceleration" of order a/Mp. This gravitational field will be negligible provided 





a k < M p 2 # 


a k > 


a k H /e k < 1 


i? ~ 3<7 k /e k 
B « sf at) 


« l/2vrGcr fc 
B « tt/GHq 


a k H /e k > 1 


u -"0 

B « 2ir 2 a k H Q 3 



TABLE I: Values of the radius i? of the Coleman-de Luccia instanton and the corresponding bounce 
action B in different limits 



that the corresponding Rindler radius (or inverse of the acceleration) is much larger than 
the radius R of the Euclidean worldsheet, which in turn is smaller than H^ 1 

a < M 2 P H . (55) 

In this regime, we can distinguish two cases. For aH Q /e <C 1, the radius of the Euclidean 
worldsheet is much smaller than the de Sitter radius, and the flat space expression (JHOJ) 
holds. In the opposite limit, aH Q /e ^> 1, we have (H^ 1 — R) <C H' 1 and 

B (wail) w 2 7i 2 aH; 3 . (56) 

The vacuum energy difference e is unimportant in this case, and the action coincides with 
that for domain wall nucleation 

Finally, the gravitational field of the brane is important when 

a > MpH Q . (57) 

In this case, the radius of the worldsheet is given by R ~ (1/2tiG(j), and 

B { W all) ^ (5g) 

In this limit, the action of the instanton is much smaller than the action of the background, 
and this is the reason why (J55J) is independent of the tension. 

The same arguments apply of course to the instantons with coincident branes, and the 
corresponding expressions for the action and radii in the different regimes are summarized 
in Table I. 



VII. THE PREFACTOR FOR CDL INSTANTONS 

The prefactor D in Eq. (|37j) is given by 

D = %-. (59) 

Here, Z and Zg are the Gaussian integrals of small fluctuations around the instanton and the 
background solutions respectively. Expanding all brane and bulk fields (which we generically 
denote by (p) around the instanton configuration (pi as (p = (pi + ^ ■ we have Se[4>] = 



Se[I] + S^ 2 ' [I, 8(f)j] + where S^ 2 ' includes the terms quadratic in 5<pj. At the one loop 
order we have 

z' = j nj&tffo) e- s<2> . 

In the functional integral, there are some directions which correspond to spacetime transla- 
tions of the instanton. The primes in the numerator of ()59)1 and in the previous equation 
indicate that the translational zero modes are excluded from the integration, and replaced 
by the corresponding spacetime volume. The latter is subsequently factored out in order to 
obtain a nucleation rate per unit time and volume. 

The degrees of freedom which live on the brane will make a contribution to the numerator 
but not to the denominator. Consider, for instance, a free bosonic field $ of mass m$ living 
on the worldsheet. Its contribution to the prefactor is 

= Z* = e- w « = Jv$ e ^~ m ^ ^ 

Here, the integral in the exponent is over the Euclidean worldsheet of the brane. If we have 
k 2 of such fields, their effect on the nucleation rate is to replace 

B k ^B k + k 2 W„, (60) 

in the naive expression for the nucleation rate, T ~ e~ Bk , a replacement which can become 
very important as we increase the number of fields. This is, in essence, the observation made 
by Feng et al. that the large number of worldsheet fields might strongly affect the nucleation 
rate. 



A. Scalars 

The Euclidean worldsheet in the Coleman-de Luccia instanton is a 3-sphere. Determi- 
nantal prefactors due to scalar fields were considered in some detail in j23[. They are given 
by 

= -Cb(-2) + {y 2 /2) ln(sinTry) / xln(smx)dx, (61) 

^ JO 

where y 2 = 1 — m\R 2 and (r is the usual Riemann Zeta function. For instance, the 
contribution of a conformally coupled scalar field can be obtained by taking m% = (3/4)R~ 2 , 
which gives 

W c = -Cb(-2) + I In 2 - -^G?(3) w 0.0638 
Hence, the effective degeneracy factor contributed by a conformal scalar field is given by 



D c w e' Wc w 0.94 < 1. (62) 

The first thing to note is that this factor is not an enhancement, but a suppression. Hence, 
the prefactor cannot simply be thought of as the exponential of an entropy. More generally, 
from Eq. (|61|1. the prefactor is independent of the expansion rate of the ambient de Sitter 
space. This has implications for the scenario proposed by FMSW, as we shall discuss in the 
next Section. 



In general, the degeneracy factor will depend on R and on the mass of the particle. For 
light minimally coupled scalars, equation (JoTj) gives 

eCa'HO 

D *~^T2 B (m«i2<l). (63) 

There can be a strong enhancement in the nucleation rate if there are very light scalar 
fields. In the limit m$ — > the factor goes to infinity. This is because a massless scalar has 
a normalizable zero mode on the sphere, corresponding to the symmetry $ — > $ + const. 
In this case, the zero mode must be treated as a collective coordinate. The nucleation rate 
is proportional to the range 5$ of the field $, because the bubbles can nucleate with any 
average value of the scalar field with equal probability |23[ 

D s {m% = 0) = lim [m $J D s (m $ )](vr J R 3 ) 1 / 2 (5<l> = e^'^ R X ' 2 8<S>. (64) 

m|->0 

Finally, for large mass the expression (jfilj) leads to 

D s « exp(vrm| J R 3 /6) {m$R > 1). (65) 

The exponent of this expression can be interpreted as a renormalization of the tension of 
the stack of branes, due to the heavy scalars living on it j2^|. Indeed, the effective potential 
for a scalar field in 2+1 dimensions in the flat space limit is proportional to m|, and the 
factor of R 3 is just due to the volume of the worldsheet. The factor D s is plotted in Fig. 
for different values of m$R. 




FIG. 5: Contribution of a scalar field to the prefactor D in the nucleation rate (|47jl. as a function 
of its mass m§ (measured in units of the inverse radius of the instanton.) At low mass, the 
enhancement is due to large phase space: bubbles can nucleate with values of the field in the range 
(|66|) . which becomes larger for smaller masses. The enhancement at large mass can be understood 
as a finite renormalization of the brane tension. 

Note that there is an enhancement both at smal and at large mass, but the two have 
very different origin. The large value of (|6l5|) for light fields can be interpreted as a phase 
space enhancement. As we shall discuss in Subsection VII-D, quantum fluctuations of fields 
living on the worldsheet of the brane are characterized by a temperature T = l/2nR. The 
corresponding fluctuations in the potential term are of order m|$ 2 ~ T 3 , which corresponds 
to a root mean squared expectation value for $ of the order 

1 



Eq. ()63|) is recovered if we insert the range ()66|) in Eq. ()64|) . The lighter the field, the 
larger is the phase space factor <5$, and the larger is the nucleation rate. For m^R 3> 1 this 
argument cannot be used, since the field does not behave as effectively massless. In that limit 
the field decouples, as it should, and its effect is felt as a renormalization of the parameters 
in the lagrangian. For a scalar field, the leading effect is to renormalize the tension of the 
brane, making it lower (See Appendix A). This causes an exponential enhancement of the 
nucleation rate. 



B. U(k) fields 

Besides scalar fields, gauge bosons and fermions may live on a brane. At low energies, 
the field content on a stack of coincident branes will be model dependent. The idea is that it 
will correspond to a gauge theory whose symmetry is enhanced when branes are coincident, 
giving rise to a large number of light species on the worldsheet. The details of the theory, 
however, will depend on whether we start from 2D branes which descend directly from the 
ten dimensions, of from higher dimensional p-branes wraped on (p-2)-cycles. They will 
also depend on the details of compactification. Rather than building a particular scenario 
from first principles, here we shall try to gain some intuition by considering a toy model 
directly in four dimensions. The degrees of freedom on the stack of k coindident branes will 
be bluntly modeled by a weakly coupled SUSY U(k) gauge theory on the 2+1 dimensional 
worldsheet. This contains a U (k) gauge field, (k 2 — 1) scalar degrees of freedom in the adjoint 
representation of SU(k), and a scalar singlet; plus the corresponding fermionic degrees of 
freedom. The action for the worldsheet fields is given by 



Ssym = ~ J Vld 3 £ 



~Tr(i^F^) + Tr(D^D^) + V(0) + ... 



(67) 



where the ellipsis indicate the terms containing fermions. Here F^ v is the field strength of the 
gauge field = A^\ a , where A a are the generators of U(k), normalized by Tr(A a Afe) = S a b/2, 
and = a A a , where a are real scalar fields (a, b = l,...,k 2 ). By analogy with the well 
known case of D-branes in 10 dimensions, we shall assume that if the coincident branes are 
flat (as in the case when there is no external four-form field), then the theory is supersym- 
metric and all degrees of freedom are massless. The scalar field is a hermitian matrix and 
can always be diagonalized by a suitable U{k) gauge transformation, = diag(</3i, 
The eigenvalues (fi, i = 1, k are then interpreted as the positions of the different parallel 
branes along an axis perpendicular to them (if the codimension of the branes were higher, 
there would be additional scalar matrices representing the positions of the branes along the 
additional orthogonal directions, but here we are interested in the case of codimension one). 
In the supersymmetric case, the potential for the scalar field vanishes, V(0) = 0. However, 
for non-flat braneSjthe displacement of the stack of branes no longer behaves as an exatly 
massless field |23|,|25|, but one which couples to a combination of the worldsheet and extrin- 
sic curvatures, as well as to the background four-form field. Also, after moduli stabilization, 
the forces amongst branes are nonvanishing, and this also contributes to the potential for 
the relative displacements of the different branes. Thus, as we shall see, there will be a 
nonvanishing potential V(0) 7^ in the physical situation of our interest. 

When the positions are not coincident, the U(k) symmetry breaks to a smaller group 
because some of the gauge bosons aquire masses m\ = (g 2 /2)(ipi — (Pj) 2 - There are always 



at least k massless vectors (corresponding to Maxwell theory on the individual branes on the 
stack) and the remaining k 2 — k have double degeneracy. For example, in the case of a single 
brane, the scalar field will represent the Goldstone mode of the broken translational sym- 
metry, associated to transverse displacement of the brane. For the case of two branes, there 
are two such scalars. One of them, ip + = ((fx + (p 2 )/\/2, will corresponds to simultaneous 
motion of both branes, and is a singlet under SU(2). The other one, (f- = (cpi — y?2)/v^2, 
will correspond to relative motion of the branes and it transforms under SU(2). When the 
two branes move apart, 0_ acquires an expectation value and two of the gauge bosons get 
a mass, breaking the symmetry U(2) — > U(l) x U(l). 

The case of interest to us is not a flat brane, but the world-wolume of the 2-brane in the 
CdL instanton, which forms a 3-sphere of radius R. In this situation, we do not expect the 
theory to be supersymmetric (in particular, corrections to the effective action will appear 
at one loop, which will be related in fact to the determinantal prefactor in the nucleation 
rate). The case of a single brane is very similar to the case of a vacuum bubble, and in that 
case we know that the transverse displacements correspond to a scalar field of negative mass 
squared [23, .25] 

m 2 + = -3R- 2 . (68) 

The origin and precise value of this mass term can be understood geometrically, since it 
leads to four normalizable zero modes which are the spherical harmonics with 1 = 1. These 
correspond to the four space-time translational modes of the instanton, which have to be 
treated as collective coordinates. This scalar field has also a single negative mode, which 
is the constant / = mode. A negative mode is precisely what is needed for an instan- 
ton to contribute to the imaginary part of the vacuum energy, and hence to contribute to 
false vacuum decay p. Integrating out the transverse displacement of the brane gives a 
determinantal prefactor of the form [23| 

D+ = °^ e u-^, (69) 

where Q = VT is the spacetime volume. The prefactor in the nucleation rate (14 7|) per unit 
time and volume is obtained after dividing by Q. The above argument neglects gravity, and 
it is a good approximation when a -C M 2 H Q . The case of strong gravity a ^> M 2 H Q is far 
more complicated, since one has to integrate over fluctuations of the gravitational field in 
the bulk, and is left for further research. 

If there are 2 coincident branes, then there are two independent transverse displacements 
corresponding to the eigenvalues tp\ and <^2- The "center of mass" displacement ip + behaves 
just like in the case of a single brane. The orthogonal combination </?_ represents the brane 
separation, and the two-brane instanton will only be relevant if this second combination 
acquires a positive mass through some mechanism, so that there is a single negative mode, 
not two, and four normalizable zero modes in total. In other words, the branes must attract 
each other. As shown in Section V, if the branes have BPS charges, they in fact tend to 
repel each other once the dilaton is stabilized, and the configuration with coincident branes 
is unstable. In this case, we do not expect that there will be any instanton representing the 
nucleation of multiple branes. 

One possible way out of this conclusion, is to assume that the charges are different from 
their BPS values, due to supersymmetry breaking effects, so that the sum Q 2 defined in 
Eq. (}4l)|) is positive. In that case, the two branes attract each other with a linear potential. 



The canonical field ip- is related to the interbrane distance d by \<f-\ = o x l 2 d. Hence, the 
interbrane potential (|45jl takes the form 



vfe>_) = gv-v>_| - + ... (to) 



y ITi 



This potential is attractive at small distances, and has a maximum at = tp r 
Q 2 a 1 1 2 1 me 1 . Classically, the branes will attract at short distances. However, there is a 
danger that they will be separated by quantum fluctuations. As we discussed in the previ- 
ous Subsection (and we will argue more at length in Subsection D) quantum fluctuations of 
fields living on the worldsheet of the brane are characterized by a temperature T = l/2nR. 
The fluctuations in the potential are of order V ~ T 3 , and these correspond to a root mean 
squared expectation value for <y?_ of the order 

<fy_ ~ a l ^ 2 /R 3 Q 2 . (71) 

The stability of the two-brane instanton requires that 5ip- <C (p m . Otherwise, unsuppressed 
quantum fluctuations take the field over the barrier and the distance between the branes 
starts growing without bound. This requires 

Q 4 » ^- (72) 

If this condition is satisfied, then ip_ is trapped near the origin, and the branes stay together. 
It can be shown that for Q 4 ^> crR~ 5 the field behaves as approximately massless in the 
range given by (|71)1. so from (J6lj) its contribution to the prefactor can be estimated as 

D - - Rl/2 ^- - (73) 

As in the case of the massive field discussed in the previous Section, this expression is 
only justified when D_ 3> 1. If Q 2 is too large, then the field will not behave as massless 
in the range (|7T|). and we expect that for Q A 3> <tR~ 5 the sole effect of the field will be to 
renormalize the coefficients of operators such as the brane tension in the classical lagrangian. 
An inconvenient feature of the linear potential (|70|l is that is non-analytic at the origin, and 
hence an explicit calculation in the limit of large slope is not straightforward. Moreover, we 
cannot write down an expression for it in terms of the matrix operator 0, but just in terms 
of its eigenvalues (fi. 

Another possibility, which is more tractable from the formal point of view, is to assume 
that there is an attractive interbrane potential which is quadratic at short distances. That 
is, as in (|T0|) with Q 2 = but with a positive coefficient in front of the second term. In 
terms of the eigenvalues <pi, which represent the displacements of the branes, we assume the 
following expression for the potential 

V(^) = m 2 +V 2 + + ^m 2 ~ Vi) 2 + (74) 

ij 

where ip + = k~ x l 2 Yli=i f% = & _1//2 Tr(0) and m 2 + = —3R~ 2 [as given in Eq. 1)68)1]. while 
m 2 _ > is a new parameter which characterizes the attractive interaction at short distances. 
In terms of the field (f>, we can write the potential as 

V(0) = m^AT^Tty) 2 + m 2 _ [Tr0 2 - /^(Tr^) 2 ] (75) 



k 2 

1 2 1 2 1 2 V^" 1 

-m + x + -m_ 2^ 



k 2 

6=2 



where in the last equality we have expanded (p = (f) a X a in the basis of generators A a , and 
we have used Ai = {2k)~ l l 2 l and TrA^ = for b = 2,...,k 2 . In the symmetric phase 
(and assuming m-R <C 1) each one of the adjoint fields <pb will contribute a determinantal 
prefactor of the form 

e Cfl'(-2) 

D -~-l72 p' (m_fl<l) (76) 

where we have used ()63j) . This shows a somewhat milder dependence in i? than in the case 
of a linear interaction between branes, given in Eq. (|73|) . but still of power law form. In 
the limit of large mass D_ m exp(7rm 3 _R 3 /6), which as discussed before amounts to a finite 
renormalization of the brane tension. 

Aside from scalars, we should also consider the contributions from gauge bosons and 
fermions. For the case of a 3-sphere, these have been studied in j^ij. For vectors of mass 
rriA, the result is 

1 fg 2 R\ , /sinh(vrm Aj R)\ f mAR 2 d , , , , , , ^ ^. , n . 

(77) 

where g is the gauge coupling, which is dimensionful in three dimensions. When the branes 
are coincident, the theory is in the symmetric phase and the gauge bosons are massless. A 
massless gauge boson gives a contribution of the form 

D A = e~ WA = gR l ' 2 e^ { - 2 \ (m A = 0) (78) 

which again behaves as a power of R. A Dirac fermion of mass m# yields the contribution 



W 9 = -logcosh(vrm J R) + vr / u 2 tanh(vra)dw + 2<^(-2, 1/2) - ~<^(0, 1/2). (79) 
4 Jo "2 

For the massless case, the R dependent terms vanish and we have 

Dq, = 2- 1 / 4 e§^ ( - 2) , (m* = 0). (80) 

which is a constant independent of the radius R. 



C. Nucleation rate 

Collecting all one-loop contributions, the prefactor in Eq. fHTj) due to the weakly coupled 
U(k) gauge theory in the unbroken phase is given by 

D = ^±{D_f-\D A ) k \D^f. (81) 
Using (JHH1), (f75 j) .(j78 j) and (J8UJ) we are led to a nucleation rate per unit volume of the form 



where is the corresponding bounce action for the nucleation of k coincident branes. Note 
that the R dependence of the prefactor is simply as a power law. Here we have used the 
form (|76|) for the scalar contribution D_, corresponding to an attractive interaction amongst 
branes which is quadratic at short distances, with a curvature of the potential characterized 
by some mass parameter m_ 2 . The prefactor in Eq. (J8*2*j) has the exponential dependence 
on k 2 which counts the number of worldsheet field degrees of freedom, while the Euclidean 
action behaves approximately linearly with k. Hence, as suggested in 0, the prefactor 
can be quite important in determining the nucleation rate. 

In the scenario proposed in it was also desirable that the enhancement in the nucleation 
rate would switch off at the present time, in order to prevent the vacuum from decaying 
further. Unfortunately, the expression (|82j) does not seem to have this property. The 
prefactor depends only on the radius R, which is itself a function of various parameters, 
such as the brane tension, the charge and the ambient expansion rate, as summarized in 
Table I. According to (|82|) . an enhancement of the nucleation rate of coincident branes will 
occur for 

R ~ Min H~\ « (g/m^f, (83) 

where we have used the results of Table 1 in the first step. Consider first the situation 
where OkM~ 2 <C H Q . Note that even in the regime when R m Hq , the dependence of 
the degeneracy factor on the corresponding dS temperature is only power law, and not 
exponential as suggested in |2j. Also, it is clear the stability of our vacuum is not guaranteed 
by the smallness of the present expansion rate. The enhancement will persist provided that 
&k/ e k <^ {g/ m -) 2 i even if the ambient de Sitter temperature vanishes. More worrisome is 
the fact that for sufficiently large k we enter the regime where o^M~ 2 3> H a . In that case, 
we have 

R ~ M 2 /a k , 

which can get smaller and smaller as we increase the number of coincident branes, eventually 
leading to a catastrophic decay rate, regardless of the value of H Q . 

The expression (JH2*j) for the nucleation rate is valid for m^R <C 1. In the opposite limit, 
m_R 3> 1, the scalars decouple, contributing a finite renormalization of the parameters in 
the action (such as the brane tension and induced Newton's constant). For completeness, 
this is discussed in Appendix A, where it is shown that the renormalization of parameters 
can have a very significant impact on the nucleation rate. 

D. Discussion: temperature of a vacuum bubble 

In Ref. 0, the prefactor in the nucleation rate for the nucleation of coincident branes 
was estimated as an entropy enhancement 

D ~ e 5 , 

where, from dimensional analysis, the entropy was estimated as S cx T 2 R 2 per field degree of 
freedom. The factor R 2 is due to the area of the bubble, and T is some effective temperature. 



2 If we assume instead a linear interaction at short distances, we should use l|73[) and the behaviour changes 
to D cx _R( 9 / 2 ) _2fe2 ) but in any case the dependence is still a power of the radius R 



Although the interpretation of the prefactor as the exponential of an entropy should not be 
taken too literally, let us try and phrase the results of the previous subsection in this intuitive 
language. 

A particle detector following a geodesic in a de Sitter space responds as if it was at rest 
in a thermal bath in flat space, at the temperature T Q = H /27r. It should be kept in 
mind, however, that the dS invariant quantum state is in fact a pure state, and hence rather 
different from a true thermal state. For instance, any two detectors in geodesic relative 
motion observe the same temperature, with a perfectly isotropic distribution. This is a 
consequence of de Sitter invariance, and is in contrast with the situation in a thermal bath 
in flat space, where moving observers detect a temperature blue-shift in the direction of 
their motion relative to the bath. 

The fields living on a nucleated brane will experience some thermal effects too. The 
bubble is embedded in a dS space characterized by a temperature T Q . The interior of the 
bubble is also a dS space characterized by a different expansion rate, with corresponding 
temperature T«. The existence of two different de Sitter spaces in contact with the brane led 
the authors of Ref. J2} to consider two different possibilities for the effective temperature of 
the fields on the brane: T = T Q and T = (TqTj) 1 / 2 . However, there is in fact no ambiguity 
in the temperature of such fields 0], which is determined as follows. 

The worldsheet of the brane is an S 3 of radius R, the Euclidean de Sitter space in 2+1 
dimensions. If interactions with bulk fields are neglected, brane fields are only sensitive to 
the geometry of the worldsheet, and do not know about the properties of the ambient space. 
In this approximation, the relevant temperature is clearly the intrinsic temperature of the 
lower dimensional de Sitter space, 

T R = l/2nR. (84) 

This conclusion remains unchanged when we include interactions with bulk fields. The 
simplest way to see this, is to consider the limiting case where gravity can be ignored and 
the nucleation takes place in a flat space. There, the ambient temperature vanishes T Q = 0, 
but the fields on the brane will feel the temperature Tr, because the nucleated brane expands 
with constant acceleration a = 1/R. An accelerating observer in the Minkowski vacuum will 
detect a Rindler temperature Tr = a/2iT, which happens to coincide with the intrinsic 
worldsheet temperature. Hence, the de Sitter vacuum in the 2+1 dimensional worldsheet is 
in equilibrium with the Minkowski vacuum in the bulk. This conclusion is quite general, and 
applies also to bubbles nucleating in de Sitter. The CdL instanton has an 0(4) symmetry 
under Euclidean rotations. This becomes an 0(3, 1) symmetry after analytic continuation 
into Lorentzian time. The quantum state after bubble nucleation is expected to inherit this 
symmetry 0, E^H^I, and the only way to achieve it is if the fields on the brane are in their 
intrinsic de Sitter vacuum, which is characterized by temperature Tr. 

Note that Tr is a relatively high temperature. The radius of the instanton is always 
smaller than H~ l (see Table 1), and therefore Tr is strictly larger than T a and Tj. Neverthe- 
less, the product k 2 Tr R 2 ~ k 2 , and hence the "entropy enhancement", is independent of 
the ambient de Sitter temperature. As shown in the previous subsections, the independence 
of the prefactor D on the ambient expansion rate is only approximate, due to the anomalous 
behaviour of light fields in de Sitter space. This introduces a dependence of the effective 
action W e ff = — log-D on the radius R of the instanton (which in turn may depend on H 
in certain regimes). This dependence, however, is quite different from the one proposed in 
0, where it was suggested that the nucleation rate of coincident branes would be enhanced 



at high and would switch off at low H Q due to the drop in ambient temperature. What 
we find instead is that, if the CdL instanton for nucleation of coincident branes really exists, 
then the corresponding degeneracy factor does not necessarily switch off. 

We shall return to a discussion of this point in the concluding Section. Before that, let us 
turn our attention to a different instanton, which may be relevant to the FMSW scenario. 



VIII. PAIR CREATION OF CRITICAL BUBBLES 

Euclidean de Sitter space is compact in all spacetime directions, and (as we just discussed) 
it behaves in some respects as a system at finite temperature. One may then ask whether 
there are instantons similar to the thermal ones in flat space. These correspond to static 
bubbles, in unstable equilibrium between expansion and collapse. 

Static instantons with 0(3) symmetry have previously been considered in a variety of 
contexts, notably for the description of false vacuum decay in the presence of a black hole (See 
e -g- HI III and references therein). The particular instanton we shall consider corresponds 



to pair creation of critical bubbles in de Sitter, and to our knowledge it does not seem to 
have received much attention in the past. This is perhaps not surprising, since its action is 
higher than the action for the maximally symmetric CdL instanton. However, if the CdL 
instanton does not exist for coincident branes, the static one may turn out to be relevant 
once the degeneracy factors are taken into account. 



A. The instanton solution 



The energy of a critical bubble is different from zero, and consequently, the metric outside 
of the bubble is no longer pure de Sitter, but Schwarzschild-de Sitter (SdS). The instanton 
is a solution of the Euclidean equations of motion, with two metrics glued together at the 
locus of the wall, which is a surface of constant r in the static chart of SdS (see Fig. EJ) • The 
metric outside is given by 

ds 2 = f a (r) dt 2 + f' 1 (r) dr 2 + rW, (85) 
where dVt 2 = d6 2 + sin 2 6d(j) 2 , and 

fo(r) = (l-™t- H 2 r 2 ) . (86) 



The metric inside is given by 

ds 2 = C 2 fi (r) dt 2 + fr 1 (r) dr 2 + rW, (87) 

where 



/, (r) = (1 - H 2 r 2 ) , 

corresponding to a de Sitter solution. The constant C is determined by the condition that on 
the brane (i.e., at r = R) the two metrics must agree, which leads to C = [f a (R) //$ (R)] 1 ^ 2 . 

The parameters M and R depend on H and Hi. Their values are determined by the 
junction conditions at the brane 



[K ab ] = -4nGa k -f ab , 



(89) 



FIG. 6: Static instanton in de Sitter space. The left figure shows the geometry induced on the plane 
r, t, while keeping angular coordinates fixed, whereas the right figure shows the geometry induced 
on the plane r, <j), keeping 9 and t fixed. The vertical direction corresponds to the coordinate r, 
common to both pictures. The cosmological horizon is at r = r+, the brane is at r = R, and r = 
is the center of the static bubble of the new phase. 



where [K a b] is the difference in the extrinsic curvature K a b = (l/2)f 1 l 2 d r g a b on the two sides 
and 'jab is the world-sheet metric. Eq. ()89|) gives rise to the junction conditions, 



[g] = -AnGa k , [g'\ = 0, (90) 

where we have introduced the new function g(r) = / 1//2 (r)/r. Using Eqs. (|8fi|) and (|88jl . we 
have 

/ 1 3GM / 1 

9o9 = + 9i9i = (91) 

Hence, using (|9T1|) . g' {R) = g[(R) = SM/inakR 1 , and then gi(R) and g (R) are easily 
obtained from Eqs. (J91j) : 

47TO" k R / 3GM\ 

^ (jR) = ^T' ^( i2 ) = ft( i2 )( 1 --R-J- ( 92 ) 

From (JHrJj) and (JHHj) we have 

^ = J?- 2 -l^- Hl 9m = J?- Hl (93) 
Inserting (|9"2|) in (f9"3*j) we finally obtain a quadratic equation for ^ (i?) = x. The solution is 

1/2 



e 3a fc 
X ~ Aa k + 16M2 + 



4a, 16M P V 2 



(94) 



where we have used if 2 — I/f = 87rGe/3 = e/3M 2 . Then the parameters M and R are given 
in terms of x by 

/T 2 = x 2 + H 2 , M = 4na k R/3x. (95) 

This concludes the construction of the instanton solution for given values of the physical 
parameters a k , H and Hi. 



The equation f (r) = has three real solutions for 27 H 2 M 2 G 2 < 1. One of them, say 
r_, is negative and the other two are positive. The two positive roots correspond to the 
black hole and cosmo logical horizons. We call them respectively r s and r + . Therefore we 
can write 

H 2 

fo (r) = (r - r_) (r - r s ) (r - r+) . (96) 

r 

Of course, in our instanton the horizon at r s is not present, since the exterior metric is 
matched to an interior metric at some r = R > r s (see Fig. For r < R the metric is just 
a ball of de Sitter in the static chart, and it is regular down to the center of symmetry at 
r = 0. In general, the size of the cosmo logical horizon is given by 

^-ts-Ot 1 )- (97) 

where we have introduced the angle 



^ = - arCtan \/27^M^- 1 ' (98) 

In the limit M — > the angle ip — > — tt/2, and H D r + — > 1. 

According to Eq (f9"2~)). at the brane we have f (R) = x 2 (R — 3GM) 2 , so the equation 
f (R) = has a double zero instead of two different roots. This means that the radius of 
the instanton will coincide with the radius of one of the horizons only in the special case 
where both horizons have the same size, r s = r + = R = 3GM. In the limit r s = r + the 

exterior metric becomes the Nariai solution ,HH, which has r = (y/3H ) \ Note that in 

the limit 3GM U/SHg) , ip — > and H Q r + — > l/y/3, as expected. 

Like in the case of instantons describing the production of black holes |3(| or monopoles 
22l | in de Sitter, the instanton presented here describes the creation of pairs of bubbles. As 



we shall see, the Euclidean solution is periodic in the time direction, so that time runs on 
a circle S 1 (See Fig. 6). The geometry at the time of nucleation is obtained by slicing the 
compact instanton through a smooth spacelike surface which cuts the S 1 factor at two places, 
say, t — and t = n. The resulting geometry contains two different bubbles separated by a 
distance comparable to the inverse expansion rate. 



B. Temperature and action 

In order to calculate the temperature of the worldsheet in the static instanton we must 
first find the time periodicity (3 . This is determined by the regularity of the Euclidean 
metric at the cosmological horizon. For r — > r + , we have 

f (r) - A 2 (l - £) , 

where 

A 2 = H 2 (r+ - r_) (r+ - r s ) = 3H 2 r 2 + - 1. (100) 
In terms of the new coordinates 



(99) 



the metric (|85|) for r — > r + reads 



ds 2 = p 2 rf0 2 + rf p 2 + r 2 dft 2 , (102) 

so it is clear that is an angle, < <p < 2tt, and t varies in the range < t < 4:7cr + /A 2 . 
Therefore the value of /3 is 

47rr + 27rr 2 
P= 3Hy + -l = r + -3GM' (103) 

The temperature of the worldsheet instanton is given by the proper time periodicity /3r = 
Jo 3 ^'° //2 (-^) ^ = /<> (-^) ^ = C fl^ 2 (R) /3. Hence, the inverse temperature is given by 

2 R - 3GM . . 

^ = 2 ™+77^GM- (1 ° 4) 

We shall also be interested in the Euclidean action, which turns out to have a rather simple 
expression in terms of r+. This is derived in Appendix B, where it is shown that the difference 
in Euclidean actions between the instanton and the background solutions is given by 

B = (1 - r 2 + Hl) . (105) 



C. Some limiting cases 

Let us start with the case of low tension branes, cr^/M 2 -C H a , H — H^. In this case the 
parameter x is large compared with H a , R ~ x _1 is small, and GMH Q <C 1. In this limit 
the angle ip in Eq. ()98j1 is close to —tt/2 and H Q r + ~ 1. We have 

„ 27r IQ-nai 2ix , 

B ^K^< hs: T; (106) 

This is just the flat space expression for the energy of a critical bubble, multiplied by 
Euclidean time periodicity of the low curvature de Sitter space in which this bubble is 
embedded. 

Next, we may consider the case of intermediate tension H Q — Hi <C cx^/M 2 -C H Q , Hi. In 
this « Hi/y/2, R « (V3x)-\ H r + « 1 - GMH a , with GMH < 1, and we have 

167T 2 (7 1. 27T , 

B "w^r h ~vw, (107) 

In this case, the difference in pressure between inside and outside of the brane is insignificant 
compared with the brane tension term, which is balanced against collapse by the cosmological 
expansion. 

Finally, in the limit of very large <Jk we find that 3GM becomes larger than R, namely, 
3GM — > 4R/3. This means that f (R) vanishes for some value = cr max , given below 
in Eq. (|118)) . so it is not sensible to consider the limit of very large but just the limit 
°k — > (Vax- As we have mentioned [see the discussion below Eq. (JHHJ)], the exterior metric 



in this limit corresponds to the Nariai solution, with r s = r + = (V3H ) . Repl acine this 
value in (J105J) we find readily 

B = mk- (108) 

It is interesting to compare this value of B with the corresponding one for the nucleation of 
black holes in the same de Sitter universe. This is described by the Nariai instanton [30], 
which has the bounce action 

B N = (109) 

Note that the difference B — = tt/3GH 2 = A^/AG, where = 4nr 2 is the area of the 
black hole horizon in the Nariai solution. Hence, the probability of nucleating black holes 
divided by the probability of nucleating brane configurations characterized by the same 
mass parameter is just the exponential of the black hole entropy, as expected from general 
considerations (In this argument, we are of course neglecting the entropy stored in the field 
degrees of freedom living on the branes, which will only show up when the determinantal 
prefactor in the nucleation rate is evaluated). 

Let us consider the value of /3 'r in the limit at — > cr max - This is a singular limit in Eq. 
(I104j) due to the simultaneous vanishing of numerator and denominator. Thus we will need 
to change to more appropriate coordinates. The fact that r s = r + does not mean, though, 
that both horizons coincide, since the coordinates r, t become inadequate in this case. Near 
this limit the metric outside takes the form ([85)1 . with 

/o(r)-^(l-£)-(l-£) 2 , (110) 

and r«ir + , plus higher orders in A. The constant A is the same parameter defined in (jlOOj) . 
but in this limit tends to zero, A 2 = V3H Q (r + — r s ). Now we define new coordinates tp and 
A by 

If r\ A 2 
008^ = 1-^1- -J, A = —t, (111) 

so that the metric becomes 

ds 2 = sin 2 ^ d\ 2 + r 2 + dip 2 + r 2 + dtt 2 . (112) 

In these coordinates the cosmological horizon is at ip = and the black hole horizon is at 
ip — 7T. Now in the limit A — » we just replace r + = (\f3H ) 

We must determine the position ipR of the brane, which is given as before by the matching 
conditions ()89j1 . where now the metric outside is ()112j) . So, on the brane, we have 

ds 2 ak = sin 2 ip R d\ 2 + r 2 + dn 2 (113) 
= f l (R)dt' 2 + R 2 dQ 2 (114) 

and the extrinsic curvature on the outside of the brane is — (l/2)d^g a b, with goo = sin 2, 
and cjqq = r 2 , i.e., K 00 = —(l/r + )g 00 cotilj,KQQ = 0. The curvature inside is as before 

K 00 = 9ood r fi /2 and K Qn = gnnf^ 2 '/r, with /• (r) = (1 - Hfr 2 ), so the Israel conditions 
give 

-^-cot^H- (fi ,2 )\ R = -4*Ga k , (115) 
f! /2 (R)/R = AnGa k . (116) 



These equations are easily solved and give 



o tt2 _ tt2 \ !/ 2 

sinfe =lera • (117) 



at = <W = 2Af p 2 J - ff, 2 ( 1 18) 

so Hi must be less than \/3H . Now regularity at the cosmo logical horizon ip ~ implies 
that < A/r + < 27r, so /3r = sin (i/)r) 2nr + . Hence, 



Thus, also in this case, the effective temperature of the field degrees of freedom living on 
the worldsheet will be of order H Q or higher. 



D. The prefactor for static instantons 

In flat space, and at finite temperature T ^> e^/afc, the instanton which is relevant for 
vacuum decay is static and spherically symmetric in the spatial directions. The fluctuations 
are periodic in Euclidean time, with periodicity (3 = 1/T. The worldsheet of the brane has 
the topology S 1 x S 2 , where the S 1 is the direction of imaginary time, and the S 2 is the 
boundary of the "critical" bubble, a closed brane in unstable equilibrium between expansion 
and collapse (this is in contrast with the zero temperature instanton, where the world-sheet 
is a 3-sphere.) The radius of the critical bubble is given by 

D ™ 

— , 

and the difference of the instanton action and the action for the background is given by 

where = (47r/3)cjfci?^ is the classical energy of the critical bubble. The one loop quantum 
correction can be written as (see e.g. 



W* = = (3{E^ - TSs) (120) 

Here, F$ denotes the free energy and E$> is the correction to the energy of the critical bubble 
due to the field $. This includes the zero point energy of $ in the presence of the bubble, 
as well as the thermal contributions. Finally, S$ is the entropy. Thus, the nucleation rate 
takes the form 

r> = De- B * ~ e ~^ +k2w ^ ~ e-? F ~ e~ E / T e s , 

where E = E^ + k 2 Eq, is the total energy, F = E^ + k 2 F§ is the total free energy and 
S = k 2 S<s> is the total entropy. Thus, for thermal instantons the determinant prefactor does 
indeed include the exponential of the entropy. This is, however, not the only role of the 
prefactor, since there is also some correction to the energy of the bubble. 



Consider, for simplicity a massless field $. At sufficiently high temperature T 3> Rp , 
the entropy behaves as S$ oc T 2 R^. From S$ = —dF^/dT, it follows that F$ « — T S$/3. 
Hence 

In this case, the prefactor clearly represents an "entropy enhancement". On the other hand, 
at lower temperatures, it is not clear whether the prefactor represents an enhancement or 
a suppression. At temperatures comparable to -R^ 1 , the vacuum energy term can be as 
important as the thermal contributions, and the logarithm of the prefactor D can have 
either sign. 

The case of the static instanton in de Sitter space is somewhat close to this low tem- 
perature situation, and without an explicit calculation it is not clear whether the prefactor 
represents an enhancement or a suppression of the nucleation rate. An interesting possibil- 
ity would be that at sufficiently high Gibbons- Hawking temperature T Q ~ H a , the thermal 
contribution may be sufficient to restore the symmetry, creating the desired attractive force 
amongst the branes. This is currently under research. 

The static instanton presented in this Section may perhaps be better suited to the sce- 
nario proposed by Feng et al. [2], than the Coleman-de Luccia instanton. Ignoring the 
degeneracy factors, the action of the static instanton is always larger than the action of 
the CdL instanton. In this sense, it seems to correspond to a subdominant decay channel. 
However, as we have discussed in previous Sections, it might well be that the CdL instanton 
for multiple brane nucleation simply does not exist because of the repulsive force amongst 
the branes. This does not exclude the possibility that in the case of the static thermal in- 
stantons the symmetry is restored at high ambient (and worldsheet) dS temperature ~ H Q . 
In this situation, the decay through nucleation of coincident branes would only be possible 
through the static instanton 3 . At low H a , the thermal contribution might not be sufficient 
to restore the symmety and stacks of branes may simply not hold together, destroying the 
possibility of further decay by coincident brane nucleation. Also, the prefactor and the 
Euclidean action have an exponential dependence on the ambient temperature ~ H a , and 
can be much suppressed at the present epoch, contributing to the stability of the present 
vacuum (in contrast with the Coleman-de Luccia case). 



E. An entropy bound 

A potentially worrying aspect of coincident brane nucleation in the CdL case is whether 
the degeneracy factor may grow without bound as we increase the number of branes 0. 
As we have seen, this will not happen for the case of the static instanton discussed in this 
section, since nucleation of coincident branes cannot involve arbitrarily large k. Indeed, 
there is a maximum value of the combined tension of the branes <C o max ~ MpH , given 
in Eq. ()118|) . beyond which the instanton simply does not exist. In this limit, the metric 
outside of the branes approaches the Nariai solution. 

From this observation, we can easily derive a bound on the maximum entropy which can 
be stored in the stack of branes. Indeed, the static instanton represents a spherical bubble 
in unstable equilibrium between undeterred expansion or collapse into a black hole. The 



3 In the weak coupling limit, we have checked that indeed the U(k) symmetry is not restored in the case of 
the CdL instanton |2J|. 



entropy can only increase when the stack of branes collapses, and hence the entropy of the 
coincident branes cannot exceed the entropy of the Nariai black hole. The stack of branes 
with tension — > <J max has the same radius as the horizon of the Nariai black hole, and 
so, the entropy of the stack of branes is bounded by one fourth of its own area, in natural 
units. A corollary is that the nucleation rate of coincident branes in the limit — > a max is 
bounded above by the nucleation rate of Nariai black holes of the same mass. 

IX. CONCLUSIONS AND DISCUSSION 

In this paper we have investigated the possibility of coincident brane nucleation by a 
four-form field, in connection with string motivated scenarios for the neutralization of the 
effective cosmological constant. 

In four dimensions, and after the moduli are stabilized, the branes repel each other at 
distances larger than the inverse mass of the moduli. At shorter distances, their interactions 
will be model dependent, but in the simplest models the branes do not attract at the classical 
level. In this situation, it is unclear whether the Coleman-de Luccia (CdL) instanton for 
nucleation of coincident branes really contributes to the semiclassical decay rate, since it 
would have too many zero modes and negative modes. 

Assuming that the CdL instanton exists for the nucleation of coincident branes (that is, 
assuming an attractive short range interaction amongst the branes in the stack), we have 
investigated the degeneracy factor accompanying the formula for the nucleation rate, due 
to the large number of worldsheet degrees of freedom. We have modeled such degrees of 
freedom by a weakly coupled SYM U(k) gauge theory, which is unbroken when the branes 
are coincident. We find that the degeneracy factor does not depend very strongly on the 
ambient de Sitter (dS) temperatures before or after the nucleation event. Rather, it depends 
only on the radius of the instanton. Hence the degeneracy factors can be quite important 
even when the ambient dS temperature is as low as it is today. This may indicate that 
nucleation of coincident branes via the CdL instanton is in fact impossible, otherwise the 
present vacuum would immediately decay. 

If the CdL instanton for coincident branes does not exist, stacks of branes may still 
nucleate through a "static" instanton which represents pair creation of critical bubbles, in 
unstable equilibrium between expansion and collapse. This is the analog of the instanton 
for thermal activation in flat space. Despite the absence of a classical attractive force, the 
branes could be held together by thermal corrections to the interbrane potential, which tend 
to favor the symmetric phase (where branes are on top of each other). The calculation of 
this thermal effective potential for the static instanton is currently under research. One may 
ask whether a similar symmetry restoration may not happen for the CdL instanton. In this 
case the calculation has been done in |24(, where it is shown that the one loop potential does 
not help restoring the symmetry. So it is conceivable that the branes may stick together for 
the static instanton but not for the CdL instanton, in which case the former would be the 
relevant decay channel. 

To conclude, we have presented some evidence that the "saltatory" relaxation scenario 
of 0] may be difficult to implement via the CdL instanton, since saltation would be hard 
to stop at present. Rather, we have speculated that it may be easier to implement through 
the static instanton. In the scenarios proposed in Ref. |2| for the saltatory relaxation of the 
cosmological constant, two different possibilities were suggested for the effective temperature 
of the worldsheet degrees of freedom, namely T\ ~ H a and T 2 ~ (H Hi) l l 2 , where Hi and 



H Q are the expansion rates before and after nucleation. We have shown that for the static 
instanton, the relevant temperature is comparable to the ambient de Sitter temperature ~ H 
before the tunneling. Hence, the nucleation rate of coincident branes would be unsupressed 
at large ambient de Sitter temperature, but exponentially suppressed at present, which is of 
course desirable. 

Clearly, many issues need to be addressed before a scenario based on coincident brane 
nucleation can be used to successfully explain the smallness of the observed cosmological 
constant. A considerable advance would be to understand why the large A e ff relaxes to 
the small A obs instead of plunging directly into deep AdS space (the latter jump involves a 
larger number of coincident branes and would be rewarded by a larger degeneracy factor). 
In Ref. an explanation was offered, based on a "uniquely weak" form of the anthropic 
principle. As explained in Section II, any relaxation mechanism requires the gap AA in 
the discretuum of A e ff not to be much larger than A a b s (otherwise it becomes a problem 
to understand why, accidentally, there happens to be an allowed vacumm so close to zero, 
at A e ff = A bs "C AA). In [2j it was proposed that AA = aA obs with a ~ 1, saturating 
the above requirement. Then the allowed A e ff would take values in the sequence (1 — 
a)A obs , A obs , (l + a)A obs , (l + 2a)A obs , .... If we start from a large A e ff, then the enhancement 
of brane nucleation for large k favours a jump to the lowest value in the above list which is 
still compatible with the existence of observers. FMSW suggested that the value (1 — a)A ofes 
may already be too small for observers to emerge, making the vacuum with the value A obs 
the favourite destination. 

Finally, one should try to embed this scenario in a cosmological context, taking into 
account the restrictions imposed by homogeneity and isotropy. If unsuppressed saltation 
happened after inflation, then we would have seen signals of it in the microwave back- 
ground. Indeed, bubbles which nucleate after thermalization are still rather small at the 
time of decoupling, and we would see different domains with different values of A e ff sepa- 
rated by fast moving stacks of branes, which would presumably cause large perturbations 
in the gravitational potential. Hence, saltation should occur during inflation, and switch 
off somewhat before the end of it. This may impose certain constraints on the space of 
parameters such as the tension and charges of the branes, or alternatively, on the ambient 
temperature below which the instanton with coincident branes simply does not exist (e.g. 
because thermal symmetry restoration is no longer effective). Also, it should be clarified 
what might be the advantages of a saltatory "neutralization" scenario over the "random- 
ization" scenarios discussed in Section II. A possible advantage is that saltatory relaxation 
operates very quickly, and hence it does not require eternal inflation to take place (as re- 
quired in the randomization scenarios). A fuller discussion of these issues is left for further 
research. 

Acknowledgements 

J.G. is grateful Roberto Emparan and Alex Vilenkin for very useful discussions. A.M. was 
supported by a postdoctoral fellowship of the Ministerio de Educacion y Cultura, Spain. The 
work by J.G. was supported by CICYT Research Projects FPA2002-3598, FPA2002-00748, 
and DURSI 2001-SGR-0061. 



APPENDIX A: HEAVY FIELDS ON THE BRANES 



The expression (|82|) for the nucleation rate is only valid for m_R <C 1. For completeness, 
here we shall discuss the limit m_i? ^> 1. In this case the heavy degrees of freedom decouple, 
and they simply contribute a finite renormalization of the parameters in front of different 
operators in the classical Lagrangian. 

For scalar fields we have 

W_ = --m 3 _R 3 + -m^R+ ... (Al) 
6 4 

As discussed in Section VII, the first term in this expansion corresponds to a finite renor- 
malization of the brane tension, multiplied by the worldvolume of the stack of branes. The 
second term correspond to a finite renormalization of the coefficient in front of the world- 
sheet Ricci scalar. This term was not present in the classical action we started with, but 
evidently it can be generated by quantum corrections. The scalar contribution (|A1J) tends 
to decrease the tension of the stack of branes. This tends to favour the nucleation rate at 
large m_R, as represented in Fig. |S| However, whether an actual enhancement really occurs 
will be model dependent, since all massive species, and not just the scalars, contribute finite 
renormalizations of the parameters in the action. In Section VII we have assumed that there 
is an attractive short range force amongst the branes, which we have modeled as a mass 
term for the scalars representing the relative positions of the branes. One may expect that 
the same mechanism which generates a potential for the scalars may generate masses also 
for their fermionic partners. From ()79jl . heavy fermions give a contribution to the effective 
potential of the form 

= + ^m%R 3 + ^m^R+ ... {m^R > 1) (A2) 
zoo 

per degree of freedom. This gives a positive renormalization of the brane tension (which 
tends to suppress the nucleation rate). In the special case where m_ = m^, the brane tension 
does not renormalize, but each pair of fields will still contribute a finite renormalization of the 
subleading term AW = (37r/8)m_i? which would suppress the nucleation rate. Generically, 
however, scalars and fermions may wind up with different masses (since supersymmetry is 
broken at some level), and the tension will be renormalized. Dividing the leading term 
in W by the worldvolume VolfS* 3 ] = 27r 2 R 3 , each scalar and fermionic degree of freedom 
contributes a brane tension renormalization of the form 

m% — m 3 . , , 

Aa « * - . (A3) 

1Z7T 

Similarly, there will be a renormalization of the induced Newton's constant Gn on the 
worldsheet, of the order AG^ 1 ~ — m_) for each pair of heavy field species. These 
changes will modify the instanton solution. For k coincident branes, the number of such 
fields grows as k 2 , and the effect of these finite renormalizations can be quite dramatic. The 
nucleation rate will take the form 

T k ~ a 2 R 2 {Ag 2 Rf' 2 e- B l cn (m_i? » 1). (A4) 

where A ~ 1 is a constant which depends on whether some (or all) of the fermionic species 
have decoupled or not. The bounce action B r k en in the exponent is calculated by using the 
renormalized values of the parameters. If Act < 0, then the renormalized tension cr& ~ 
ka + ik 2 — 1)A<7 sharply decreases for large k, leading to unsuppressed tunneling rate. On 
the contrary, for Act > the nucleation of coincident branes is strongly suppressed. 



APPENDIX B: EUCLIDEAN ACTION FOR THE STATIC INSTANTON 



The action is given by 



S E (I) = a J d 3 ^ + J d'x^g (p v - . (Bl) 

On shell, the scalar curvature is given by, 

TZy/g = 32irGp Vy /g + 24nGa J d 3 ^S (4) (x - x (£)) , (B2) 
and hence the instanton action is given by 

S E (I) = ~J d 3 ^- J d 4 xp v ^. (B3) 



The first integral in ()B3j) is just the volume of a two-sphere of radius R times (3r. The 
second integral in (|B3)1 splits into two contributions from the two regions, 



pR rr+ 

Pi / CdtdrA-nr 2 + p Q dtdrA-nr 2 (B4) 
Jo Jr 

= Pl C(3^R 3 + Po (3^k (r% - R 3 ) (B5) 

So the instanton action is 

S E (I) = -2nR 2 a f y* (R) (3 - R^j^P ~ K " ^) ^ (B6) 
After some algebra S E (I) can be written in the simple form 

7TT 2 

Se(I) = (B7) 

The exponent B which gives the probability for brane nucleation is the difference in Eu- 
clidean actions between instanton and background. The action of the background is just 
Se = —^/GH 2 , so the difference in Euclidean actions leads to Eq. (jl()5|) . 
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